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ABSTRACT 



X-ray colour-colour (CD) and colour-intensity (HID) diagrams are powerful tools that allow the investigation of spectral variability 
without the assumption of any spectral model. These diagrams have been used extensively in low-mass X-ray binaries and black-hole 
candidates but very few applications are found for high-mass X-ray binaries. We have investigated the spectral and timing variability 
of four accreting X-ray pulsars with Be-type companions during major X-ray outbursts. The aim is to define source states based on 
the properties (noise components) of the aperiodic variability in correlation with the position in the colour-colour diagram. Different 
spectral states were defined according to the value of the X-ray colours and flux. Transient Be/X-ray binaries exhibit two branches 
in their colour-colour and colour-intensity diagrams: the horizontal branch corresponds to a low-intensity state and shows the larger 
fractional rms, similar to the the island state in atolls and horizontal branch in Z sources; the diagonal branch corresponds to a high- 
intensity state, in which the source spends about 75% of the total duration of the outburst. Despite the complexity of the power spectra 
due to the peaks of the pulse period and its harmonics, the aperiodic variability of Be/X-ray binaries can be described with a relatively 
low number of Lorentzian components. Some of these components can be associated with the same type of noise as that seen in 
low-mass X-ray binaries, although the characteristic frequencies are about one order of magnitude lower. The analysis of the CD/HID 
and power spectra results in two difi'erent types of Be/X. While in 4U 0115+63, KS 1947+300 and EXO 2030+375 the hard colour 
decreases as the count rate decreases, in V0332+53 it increases. The pattern traced by V0332+53 then results in a Z shaped track, 
similar to the low-mass Z sources, without the flaring branch. In contrast, the horizontal branch in 4U 0115+63, KS 1947+300 and 
EXO 2030+375 corresponds to a low/soft state, not seen in other types of X-ray binaries. The noise at very low frequencies follows 
a power law in V0332+53 (like in LMXB Z) and it is flat-topped in 4U 01 15+63, KS 1947+300 and EXO 2030+375 (like in LMXB 
atoll). V0332+53 shows a noise component coupled with the periodic variability that it is not seen in any of the other three sources. 



Key words, stars: individual: 4U 0115+63, V 0332+53, KS 1947+300, EXO 2030+375 
binaries close -stars: emission line. Be 



X-rays: binaries - stars: neutron - stars: 



1. Introduction 

, The colour-colour (CD) and hardness-intensity (HID) diagrams 
are very useful tools to investigate the spectral variations of an 
X-ray source. They are model- and instrument-independent and 
reflect the intrinsic properties of the system. The first application 
of this type of plots was the separation of different types of galac- 
tic X-ray sources and the classification of X-ray populations. The 
spectral hardness was defined by the ratio of counting rates in 
each of the proportional counters of flie MIT OSO-7 experiment 
('Marker tetai]|l977l: lRappaportet"anil977l) . IWhite & Marshalll 
( 1984) made use of the X-ray colour-colour diagrams, obtained 
from HEAO 1 observations, to study the spectral states among 
the different types of X-ray binaries, including black-hole candi- 
dates. 

CDs became widely used during the EXOSAT era, which 
also marks the introduction of HIDs (Priedhorskv et al. 1986; 
[van der Klis et al.l Il987al; iHasingerl il987). where "hardness" 
refers to the ratio of the count rate in two energy bands. Before 
EXOSAT, CDs contained sources of different origin. The nov- 
elty of the EXOSAT CDs and HIDs was that they were applied 
to individual sources. CD/HIDs became a useful way to inves- 
tigate the rapid aperiodic variability in correlation with spectral 
states, as defined by the spectral hardness, of individual low- 
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mass X-ray binarie s and black-hole candidates dSchulz et al.l 
ll989l;lHasinger & van der Klis 1989) . 

With the advent of RXTE flie use of CDs and HIDs have 
become part of the standard data analysis. RXTE has provided 
an unprecedentedly large database of X-ray observations with 
exceptionally good timing and moderate spectral resolution. 
Together with the functional form of the variability components 
("noise"), CD and HID have become an essential tool in devel- 
oping a whole new phenomenology on the spectral and tim- 
ing properties of X-ray binaries by introducing the notion of 
source states. A state is defined by the appearance of a spec- 
tral (i.e. power-law, blackbody) or variability component (i.e. 
Lorentzian) associated with a particular and well-defined posi- 
tion of the source in the CD/HID. A transition between states 
is assumed when the relative strength of the spectral or vari- 
ability component s varies and the sou rce motion in the CD/HID 
changes direction dvan der KIisll2006l) . 

While there are numerous references in the literature on 
the application of C D/HID analysis on low-mass X-ray binaries 
dvan der KIisll2006|, and reference s therein) and black-hole sys- 
tems (see e.g. iBelloni et alj|2005l) . very litt le work of this type 
has been done on high-mass X-ray binaries. iBelloni & Hasingen 
(1990) performed an aperiodic variability study of 12 high- 
mass X-ray binaries from EXOSAT archive data. The sample 
included EXO 2030+375 and V0332+53. However, the few 
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Table 1. Optical and X-ray information of the systems and the outbursts analysed in this work. 



Source 


Spectral 


i^spin 


Porb 


e 


Distance 


Outburst 


ix.pcak 


On-source 


name 


type 


(s) 


(days) 




(kpc) 


duration (d) 


3-30 keV (erg s"') 


time (ks) 


4U 0115+63 


B0.2V 


3.6 


24.3 


0.34 


8.1 


55 


1.2x10'** 


91 


KS 1947+300 


BOV 


18.7 


40.4 


0.03 


10 


165 


6.2 X 10" 


132 


EXO 2030+375 


B0.5III-V 


41.8 


46.0 


0.41 


7.1 


155 


1.4 xlO'* 


166 


V 0332+53 


08-9V 


4.4 


34.2 


0.30 


7 


105 


3.2 X lO'** 


198 



observations analysed of V0332+53 did not correspond to a 
major X-ray outburst, while the data of EXO 2030+375 cov- 
ered the decay of the 1985 outburst only partially. The other 
10 sources were supergiant X-ray binaries in bright states. The 
vast majority of the rapid aperiodic variability studies in accret- 
ing X-ray pulsars have concentrated on the detection of QPOs 
dAngeHni et al."l989'; 'Jernigan et al.||2000t |Qu et al]|2005|); see 
also Shirakawa & Lai (2002, and references therein) 

Massive X-ray binaries are classified according to the lu- 
minosity class of the optical component into Be/X-ray binaries 
(dwarf or subgiant) and supergiant (luminosity class I-II) X-ray 
binaries. The former tend to be transient systems while the latter 
are persistent sources. Be stars are non-supergiant fast-rotating 
B-type and luminosity class III-V stars which at some point 
of their lives have shown spectral lines in emission. In the in- 
frared they are brighter than their non-emitting counterparts of 
the same spectral type. The line emission and infrared excess 
originate in extended circumstellar envelopes of ionized gas sur- 
rounding the equator of the B star When the Be star takes part 
on a binary system where the companion is a neutron star, then 
the system is referred to as a Be/X-ray binary. 

Be/X-ray binaries are variable on time scales from seconds to 
years. The fastest variability is found in the X-ray band. Virtually 
all the Be/X-ray binaries with identified optical counterparts are 
X-ray pulsars. Pulse periods cover the range 1-10^ s. On longer 
time scales (months to years) the variability is also apparent in 
the optical and infrared bands and it is attributed to structural 
changes of the circumstellar disc. The long-term X-ray variabil- 
ity of the transient Be/X-ray binaries is characterised by two type 
of outbursting activity 

- Type I outbursts. These are regular and (quasi)periodic out- 
bursts, normally peaking at or close to periastron passage of 
the neutron star. They are short-lived, i.e., tend to cover a rel- 
atively small fraction of the orbital period (typically 0.2-0.3 
Poib)- The X-ray flux increases by about one or two order 
of magnitude with respect to the pre-outburst state, reaching 
peak luminosities Lj < lO-'^ erg s"'. 

- Type II outbursts represent major increases of the X-ray flux, 

10^-10"* times that at quiescence. They reach the Eddington 

luminosity for a neutron star (~ lO"*^ erg s"') and become the 

brightest objects of the X-ray sky. They do not show any pre- 

feiTed orbital phase and last for a large fraction of an orbital 

period or even for several orbital periods. The formation of 

J II 1 

an accretion disc during Type II outbursts (Kriss et al .lll983t 

Motchetal 1991; Hayasaki & Ozakaki 2004; Wils on et alj 
2008) may occur The discovery of quasi-periodic oscilla- 
tions i n some systems ( | Angelini et alJl989tlTakeshima et all 
119941 ; iFmgeretalJI 19961) would support this scenario. The 
presence of an accretion disk also helps explain the large and 
steady spin-up rates seen during the giant outbursts, which 
are difficult to account for by means of direct accretion. 

In this work we present a systematic study of the aperiodic 
variability of four Be/X-ray binaries during type II outbursts. 



First, we give a summary of the relevant information of the 
sources. Sect. [3] summarises the observations and the character- 
istics of the instruments used in the analysis. Sect.|4]refers to the 
reduction of the data and describes how the CD/HID and power 
spectra were obtained. A detailed description of the results is 
given in Sect. |5] This section is divided into three subsections 
describing the outburst profiles, the colour analysis and the ape- 
riodic variability. Our interpretation of the results and the impli- 
cations of this work are discussed in Sect. |6] Final conclusions 
are drawn in Sect. [T] 



2. The sources 

2A. 4U 0115+63 

4U 0115+63 is one of the most active and best studied Be/X- 
ray transients as it was one of the first Be/X-ray binaries to be 
discovered. The oldest available X-ray observations date back to 
August 1969 when the Vela 5B satellite detect ed the source a s 
three small outbursts separated by 180 days (Whitlock'1989^ 
Since then about 15 outbursts have been reported (ReigetaTj 
2007). During the RXTE's life time 3 major outbursts have been 
detected: in 1999, 2000 and 2004. Here we analyse observa- 
tions from the 2004 event. The neutron star orbits a B0.2Ve 
in a moderately eccentric orbit with Pp b - 2 4.3 d, e = 0.34 
and fl I sin/ = 140.1 It-s (Rappapo rt et al.|[l978l) and rotates with 



pin 



3.6 s, as infeiTed from the pulsed X-ray em ission. The 
estimated to be ~8 kpc jReig et al.l 



distance to the source is 
2007). 

The X-ray spectrum shows up to five cyclotron resonance 
scattering features (fundamental and four harmonics) and con- 
stitutes the ac creting neutron sta r with the highest number of cy- 
clotron lines (J Heindl et al.ll2004l) . The optical and infrared emis- 
sion is characterised by cyclic changes with a period of ~ 5 
years. Another peculiarity is that X-ray outbursts in 4U 01 15+63 
appear to come in pairs, i.e., two in ever y cycle. However , 
sometimes the second outburst is missing (iReig et al.l l2007h . 
mHZ QPOs were detected during the 1999 and 2004 outbursts 
(Heindletal.,,1999 ; Co burn et al..l2004, ). 



2.2. KS 1947+300 

The Be/X-ray binary KS 1947+300 was first detected in the 
X-ray band on June 8, 1989 by the TTM coded mask X-ray 
spectrograph on the Kvant module of the Mir space station 
(iBorozdin et alJ 1 19901) with a peak flux of 70 mCrab. About 
35 days later the flux was at the limit of sensitivity of the in- 
strument with a 3cr upper limit of 10 mCrab. X-ray pulsations 
with a pulse period of 18.7 s were found in the X-ray flux of 
GRO J 1948 +3 2 in April 1994 by the a ll-sky monitor BATSE 
aboard CGRO dChakrab artv et alj |1995h . The source reached a 
peak pulse flux of 50 mCrab in the energy range 20-75 keV. 
About 25 days later the flux decayed bellow the detection thresh- 
old of BATSE. Given the very large error box of GRO J 1948 +32 
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Table 2. Observations and definition of source states. 



Source 


Branch" 


Time range 


Proposal 


On-source 


<5C>'' 


< HC >" 


Flux^ 


Luminosity" 


state 




MJD 


ID 


time (s) 










4U 0115+63 


Risel 


DB 


5S254.0676-53256.4660 


90089 


3872 


2.08±0.05 


0.71+0.03 


0.95 


7.5 


Rise2 


DB 


53260. 1495-5S265.0921 


90089 


3536 


1.91+0.05 


0.81+0.04 


1.5 


11.8 


Peak 


DB 


53267.1686-53272.9721 


90089 


29744 


1.87±0.01 


0.89+0.01 


1.4 


11.0 


Decay 1 


DB 


53274.8282-53280.7687 


90089 


9712 


1.95±0.02 


0.85+0.02 


1.0 


7.9 


Decay2 


DB 


53278.6230-53280.7687 


90089 


6128 


2.03±0.04 


0.70+0.06 


0.8 


6.0 


Decays 


DB 


53282.5715-53282.6673 


90014 


5536 


2.09±0.01 


0.62+0.01 


0.6 


5.0 


Decay4 


HB 


53284.6063-53304.5730 


90014 


32500 


1.99±0.04 


0.57+0.03 


0.11 


0.9 


KS 1947+300 


RiseO 


HB 


51874.1502-51878.3595 


50425 


16720 


1.57±0.07 


1.24+0.02 


0.01 


0.12 


Risel 


DB 


51883.8286-51892.2963 


50425 


14060 


1.68±0.04 


1.13+0.05 


0.09 


1.0 


Rise2 


DB 


51896.1987-51910.9382 


50425 


17744 


1.55+0.03 


1.14+0.02 


0.25 


2.3 


Rises 


DB 


51914.1750-51932.3484 


50425 


10080 


1.49+0.03 


1.16+0.03 


0.36 


4.3 


Peak 


DB 


51941.6795-51968.5845 


50425 


23344 


1.40+0.02 


1.20+0.02 


0.53 


6.3 


Decay 1 


DB 


51970.2073-51990.1352 


60402 


15536 


1.51+0.03 


1.19+0.02 


0.39 


4.7 


Decay2 


DB 


51992.1719-52010.8263 


60402 


9888 


1.59+0.02 


1.16+0.02 


0.22 


2.6 


Decays 


DB 


52014.3856-52050.1873 


60402 


9984 


1.65+0.04 


1.05+0.07 


0.11 


1.3 


Decay4 


HB 


52052.9510-52078.1324 


60402 


14720 


1.32+0.08 


1.11+0.04 


0.008 


0.1 


EXO 2030+375 


Risel 


DB 


53908.6156-53924.7247 


91089 & 92067 


43200 


1.39+0.07 


1.09+0.02 


1.4 


8.4 


Rise2 


DB 


53926.6282-53943.2576 


91089 


26512 


1.26+0.02 


1.12+0.01 


2.1 


12.7 


Peak 


DB 


53944.0532-53975.4312 


91089 


35536 


1.23+0.02 


1.13+0.02 


2.4 


14.5 


Decay 1 


DB 


53976.1186-53998.6782 


91089 


25168 


1.32+0.05 


1.09+0.03 


1.6 


9.7 


Decay2 


DB 


53999.5902-54015.4534 


92067 


8656 


1.47+0.05 


0.99+0.02 


0.57 


3.4 


Decays 


DB 


54016.5928-54030.9436 


92067 


10912 


1.43+0.03 


1.02+0.04 


0.11 


0.69 


Decay4 


HB 


54031.7791-54047.9337 


92067 


16144 


1.41+0.03 


0.96+0.03 


0.06 


0.37 


V0332+53 


Risel 


DB 


53340.2869-53346.7436 


90089 


48752 


2.18+0.02 


0.89+0.03 


2.3 


13.5 


Rise2 


DB 


53352.7706-53361.1469 


90089 


32720 


2.06+0.03 


0.73+0.02 


5.0 


29.3 


Peak 


DB 


53364.9113-53369.5965 


90089 & 90427 


13040 


1.95+0.02 


0.66+0.01 


5.5 


32.2 


Decay 1 


DB 


53374.0267-53381.6610 


90427 & 90014 


23504 


2.04+0.02 


0.74+0.01 


4.2 


24.6 


Decay2 


DB 


53385.8724-53395.4413 


90014 


38672 


2.11+0.03 


0.83+0.03 


2.6 


15.2 


Decays 


DB 


53398.4560-53414.0726 


90014 


21296 


2.20+0.04 


0.94+0.03 


1.5 


8.8 


Decay4 


HB 


53416.0726-53424.4139 


90014 


20160 


2.05+0.15 


1.13+0.05 


0.29 


1.7 



a: HB: Low-intensity (Horizontal) Branch; DB: high-intensity (Diagonal) branch 

b: Normalised to the Crab. Average Crab values are < SC >crab= 0.34 and < HC >crab= 

c: xlO"* erg cm"^ s"' in the 3-30 keV band 



0.47 



d: xlO" erg s"' in the 3-30 keV band 



no connection with this source and KS 1947+300 was realised. 
The source remained i n quiescence until the beginning of 2001 
JGal loway et alj l2004t) . In this occasion the source reached a 
maximum flux of 120mCrab, the largest in its history. This 
outburst is the object of the present work. Prior to this ma- 
jor outburst, weak emission (~ 20 mCrab) and X-ray pulsa- 
tions (18.7579+0.0005 s) from KS 1947+300 had been detected 
in 2 000 November ( Levine & Corbet! 120001: ISwank & Morgan 
12000 ). The coincidence in the value of the pulse period and 
the location of KS 1947+300 inside the GRO J1948+32 error 
circle suggest that the two sources are in fact the same object. 
Pulse timing analysis allowed to solve for the orbital parameters 
resulting in orbital period Porb = 40.415 ± 0.010 d, eccentric- 
ity e = 0.033 + 0.013 and projected radius Oj^sin/ - 137 + 
3 It-s dOallowav et al.l 120041) . A broad-band spectral study of 
KS 1947+300 using BeppoSAX observations (i Naik et alj|2006l) 
showed that the energy spectrum in the 0. 1-100 keV energy band 
has three components: a Comptonized component, a ~0.6 keV 
blackbody component, and a narrow, weak iron emission line 
at 6.7 keV with a low column density of material in the line of 



sight. Although the optical counterpart to KS 19 47+300 was cor- 
rectly identified soon afte r the Kvant detection dGoranskii et al.l 
1991:'Grankin et al.'l99l') its identification as a Be/X-ray binary 
was suggested by Negueruela et al. (2003). KS 1947+300 is as- 
sociated with a moderately reddened V=14.2 BOVe star located 
at ~ lOkpc. 



2.3. EXO 2030+375 

EXO 2030+375 w as discovered by EXOSAT in 1985 
dParmar et alj Il989l) . Since then it has been extensively ob- 
served in the X/y-ray band by various missions: EXOSAT 
(Reynolds^ al. 1993), ROSAT ( Mavroma takis 199-4 
RXTE (JReig&Cod [1998L [l999l). CG RO (TStoUber et all 
[T999I). INTEGRAL dCamero-Arranz et all l2005h and SWIFT 



dKlochkov et"an l2007l) . EXO 2030+375 is die prototype of 
X-ray variability in Be/X-ray binaries. It shows a regular pattern 
of type I outbursts (Lx ~ 10^^ erg s"') during each periastron 
passage dWilson et al .1 12002). In addition, EXO 2030+375 has 
shown two major outbursts (Lx ^ 10^** erg s"'): the first one 
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Fig. 1. The outburst profiles as seen by ASM RXTE. Vertical lines mark the time of the RXTE PCA observations. Different colours 
distinguish the spectral regions. Bin size is 1 day. 1 Crab corresponds to ~75 c/s. See the electronic edition of the Journal for a 
colour version of this figure. 



in 1985 that led to its discovery as an X-ray transient and the 
second one in 2006. In this work we analyse RXTE observations 
of the latest event. 

EXO 2030+375 contains a rotating neutron star (f spin - 
41.7 s) orbiting a Be prima-y jMotch & Janot-Pachecol [l987t 
ICoe et all Il988t iReig et al] 1 19981) in a moderately eccentric 
(e - 0.4) wide (f orb = 46 d) orbit and is located at about 7.1 
kpc dWilson et alj |2002'). Quasi-periodic oscillations with a fre- 
quency of 0.2 Hz were detected from this source when the lu- 
minosity was close t o its maximum during the 1985 outburst 
dAngelini et al.|[l989h . In common with other high-mass X-ray 
binaries, the continuum spectral shape in the range 1-30 keV 
can be represented by a power law with an exponential cutoff 
dRevnolds et al. |[l993h . A blackbody component has been re- 
erg s~') lumi- 
10^*^ erg s-i) 



ported to give good fits at very high (Lx ~ 10^** 
nosity ('Sun et al."l994') and relatively low (L^ ^ 
luminosity (Reig & Coe. .l999.) . 



2.4. V0332+53 

A similar an alysis to the one presented here has already been 
published in iReig et all (|2006|) . For the sake of clarity and in 
order to apply the same reduction procedure t o all sources, we 
have re-analysed the observations in lReig et al.l (2006). Note also 
that the present work includes data during the rise of the outburst 
that were not previously analysed and that the frequency range 
of the power spe ctra is two orders of magnitude larger than in 
IReig et al.1 (l2006l) . 



V0332+53 is a hard X-ray transient that spends most of its 
life in a quiescent X-ray state. After more than 15 years of qui- 
escence the Be/X-ray binary V0332+53 underwent a giant out- 
burst in December 2004. The outburst was predicted about 9 
months before its occurrence from the brightening of the opti- 
cal companion ( Goranskii & Barsukova 2004). Prior to the event 
reported here, V0332+53 had went into outburst is three other 
occasions. The 1973 outburst led to its discovery as a bright 
X-ray source by the Vela 5B satellite dTerrell & Priedhorskvl 
1984). In 1983 it reappeared in the form of three small outbursts 
(Tanaka 1983). EXOSAT observations of this activity period re- 
sulted in the discovery of X-ray pulsations with Pspin = 4.4 s 
and the determina tion of the orbital parameters, Porb = 34.25 
d and e = 0.31 dStellaet al.lll985l) . In 1989 Ginga detected 
V0332+53 again and allowed the discover y of a cycl otron reso- 
nant scattering feat ure at 28.5 keV (Makishima et al.i 199 0) and 
QPOs at 0.051 Hz dTakeshima et al.lll994 . The optical counter- 
part to V0332+53 is an 08-9Ve star at a distance of ~ 7 kpc, 
show ing Hg in emission and strong and variable infrared emis- 
sion (Bernacca et al. '19844 ICorbetet al.1 1 19861; ICoe et alJ[T987l; 
[Negueruela et al. 1999i) . 

3. Observations 

We have analysed data obtained by the Proportional Counter 
Array (PCA) and the All Sky Monitor (ASM) onboard the Rossi 
X-ray Timing Explorer (RXTE). The ASM data consist of daily 
flux averages in the energy range 1.3-12.1 keV. The PCA cov- 
ers the lower part of the energy range 2-60 keV, and consists 
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Fig. 2. PCA light curves and hardness ratio evolution. Each point represents a 5 12-s integration. Symbols are as follow: pre-outburst 
(open black circles, only in KS 1947+300), risel (red circles), rise2 (orange squares), rise3 (pink crosses, only in KS 1947+300), 
peak (black triangle-up), decayl (violet stars), decay2 (brown diamonds), decay3 (green triangle-down) and decay4 (blue pluses). 
See the electronic edition of the Journal for a colour version of this figure. 



of five identical coaligned gas-filled proportional units. Due to shall refer to each of these data intervals as "pointings". Data 

RXTE's low-Earth orbit, the data consist of a number of contigu- taken during satellite slews, passage through the South Atlantic 

ous data intervals (typically 1 hr long) interspersed with obser- Anomaly and Earth occultation were removed, 
vational gaps produced by Earth occultations of the source and During the RXTE life time the response of the detectors var- 

passages of the satellite through the South Atlantic Anomaly. We ied due to ageing. Also, gain changes are applied occasionally. 
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Fig.3. Colour-colour diagrams of the four sources investigated in this work. Each point represents a 512-s integration. SC=7-10 
keV/4-7 keV, HC= 15-30 kev/10- 15 keV. Symbols are as follows: riseO (open black circles, only in KS 1947+300), risel (red circles), 
rise2 (orange squares), rise3 (pink crosses, only in KS 1947+300), peak (black triangle-up), decayl (violet stars), decay2 (brown 
diamonds), decay3 (green triangle-down) and decay4 (blue pluses). See the electronic edition of the Journal for a colour version of 
this figure. 



making the channel boundaries for a given energy range change 
with time, and also slightly affecting the effective areas of the de- 
tectors. Each gain change is the start of a new "gain epoch". To 
minimize this effect, data from the same "epoch" was used. All 
the data analysed in this work correspond to "epoch" 5. In ad- 
dition, the source X-ray colours were normalised to those of the 
Crab that are closest in time. Table |2] gives relevant information 
about the observations. 



4. Data reduction and analysis 

In this section we describe the different tools (CD/HID and 
power spectra) used in the data analysis. ASM light curves were 
obtained to study the outburst profiles. PCA data can be col- 
lected and telemetered to the ground in many different ways de- 
pending on the intensity of the source and the spectral and tim- 
ing resolution desired. In this work we used Standard! to obtain 
colour-colour and hardness-intensity diagrams and time-average 
energy spectra. In addition, various high-time resolution modes 
(goodjcenon, EventJ25fis_64MJ)_ls and SB_62fis_0_49) were 
used to extract the light curves with which the power spectra 
were produced. 



4.1. Outburst profiles 

The outburst profiles were obtained from the 1-day binned ASM 
light curves and are shown in Figure [1] Table [T] gives the dura- 
tion and peak intensity of each individual outburst together with 
some relevant information about the systems. The time of the 
PCA observations are marked in Fig. [1] with different colours 
indicating different spectral regions (see the electronic edition 
of the Journal for a colour version of this figure). Each vertical 
mark corresponds to a pointing. In the energy range of the ASM, 
1 Crab corresponds to ~75 c s '. In general, the PCA sampling 
of the outbursts is excellent. The only part of the outburst not 
covered by the PCA data is the 7-10 first days, which can be at- 
tributed to the transient nature of the sources (it takes some time 
to issue an alert warning and reschedule the observations). The 
only source with available data covering the entire outburst is 
KS 1947+300. 



4.2. Colour-colour arid hardness-intensity diagrams 

Colour-colour (CD) and hardness-intensity (HID) diagrams are 
powerful tools that allow the investigation of spectral variability 
without the assumption of any spectral model. In order to study 
the evolution of the spectral shape as the outburst evolves, the 
RXTE light curve was divided into different parts depending on 
count rate and whether the outburst was on the rise, peak or de- 
cay. Given the good sampling of the PCA observations, these 
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Fig. 4. Hardness-intensity diagrams. Each point represents a 5 12-s integration. Intensity corresponds to the energy range 4-30 keV. 
Symbols are as follows: riseO (open black circles, only in KS 1947+300), risel (red circles), rise2 (orange squares), rise3 (pink 
crosses, only in KS 1947+300), peak (black triangle-up), decayl (violet stars), decay2 (brown diamonds), decay3 (green triangle- 
down) and decay4 (blue pluses). See the electronic edition of the Journal for a colour version of this figure . 



three phases were further subdivided into shorter intervals. Thus 
each interval defines a spectral region characterised by count 
rate, soft and hard colours. Table |2] gives relevant information 
about each region and source. 

The CD and HID were obtained from the background- 
subtracted Standard 2 PCA light curves in the energy ranges 
ci =4-7keV,C2 =7-10keV,C3 = 10-15keVandc4 = 15-30 
keV. The soft colour (SC) was defined as the ratio C2/C1 and the 
hard colour (HC) as the ratio C4/C3. These ratios are expected 
to be insensitive to interstellar absorption effects. The hydrogen 
column densities to the systems, obtained from model fits to the 
X-ray spectra, are in the range 0.6 - 3 x 10^^ cm"^. The evolution 
with time of the PCU2 count rate and the soft and hard colours 
is shown in Fig.|2l 

The CD was constructed by plotting the hard colour as func- 
tion of the soft colour (Fig. O and the HID by plotting the count 
rate in the 4-30 keV band as a function of one of the colours 
(Fig.m. To recover the true values of the colours of the sources, 
the soft colour should be multiplied by 0.34 and the hard color 
by 0.47 (quoted values are averages of the Crab colours during 
the observations). The variation of the Crab colours computed as 
the root-mean-square, i.e. the standard deviation over the mean 
colour, was, on average, 2.1% and 1.6% for the soft and hard 
colours, respectively. 

To reduce the uncertainty in the individual 16-s points (de- 
fault resolution of the Standard 2 mode) the light curves were 
rebinned to 5 12-s bins. Also, data for which the resulting rela- 
tive errors were larger than 5% were excluded. When the count 
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rate per PCU was lower than 10 c s \ relative errors up to 10% 
were allowed. 
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Table 3. Power spectral parameters (centroid frequency, PWHM and rms) for the Lt,, Li and L„ components. Errors are 90% 
confidence level. The centroid frequency of the zero-centred Lorentzian was fixed during the fit. HB stands for horizontal branch. 
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*: including components in Table |4] and |5] 

**: affected by the mHz QPO. See Fig.[5]and Table|5] See also Cobum et al. (2004) 



Table 4. Power spectral parameters (centroid frequency, FWHM 
and rms) for the L„' component in EXO 2030+375. Errors are 
90% confidence level. 
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4.3. Power spectral analysis 

Power spectra were created by performing a fast Fourier trans- 
form on the light curves according to the following procedure: 
one light curve in the energy range 2-20 keV (channels 0^9) 
was extracted for each spectral region with a time resolution of 



2"^ seconds. Each light curve was divided into 512-s segments 
and an EFT was calculated for each segment. The final power 
spectra resulted after averaging all the individual power spectra 
and rebinning logarithmically in frequency. All power spectra 
are presented in the rms normalisation, where the power inte- 
grated over a certain frequency interval equals the squared rms 
of the source (Belloni & Hasinger 1990; Miyamoto et al. 1991). 

The plots of the power spectra are shown in the v x Py rep- 
resentation, where each power is multiplied by the correspond- 
ing frequency. This representation helps visualize at which fre- 
quency the contribu tion to the total rms variability is maximum 
dBelloni et alj 12002 ). Note that the frequency at which maxi- 
mum power is attained in the v x Py representation does not 
equal the centroid frequency of the Lorentzian, vq, but Vmax = 
{vl + {FWHMIlf)^^^, where FWHM is the Lorentzian full- 
width at half maximum. 

The statistical Poisson noise {Pm - 2/Rt, where Rj is the 
source count rate) was modified by dead time effects. Dead-time 
effects are expected to affect the power spectra at high frequen- 
cies, especially during the high-flux states, where the count rate 
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Table 5. Power spectral parameters (centroid frequency, FWHM and rms) for the L„iqpo, Lif and L, components. Errors are 90% 
confidence level. The centroid frequency of the zero-centred Lorentzian was fixed during the fit. HB stands for horizontal branch. 
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is over ~ 1000 c s ' . One way of checking whether the dead-time 
effects are important is by measuring the power at high frequen- 
cies in a Leahy normalised powe r spectrum. The expected power 
density should be at a level of 2 ( jLeahv et al.ll 19831 : Ivan der Klisl 
[1989; Jernigan et al. 2000). As an example, the power density of 
EXO 2030-1-375 during the peak of the outburst ("peak"), calcu- 
lated for V > 100 Hz, is 1 .935, while that obtained during the end 
of the outburst ("decay3") is 1 .995. The correction of the Poisson 
noise Pn{v) for dea d-time effects was calculated as described in 
iNowaketal.] ( 119991) . 

The study of the rapid aperiodic variability in X-ray pulsars 
is hampered by the peaks due to the periodic pulsations. The spin 
period and its harmonics show up as narrow peaks in the power 
spectrum that distort the continuum. The width of these peaks 
depends on the frequency resolution, which in turn, depends on 
the length of the segments (Av = l/T, where T is the total dura- 
tion of the segment). On one hand, a high frequency resolution 
is required in order to have well defined peaks. On the other 
hand, higher frequency resolution implies less number of power 
spectra to average and worse signal-to-noise, especially at higher 
frequencies. As mentioned above we found a good compromise 
with r = 512 seconds. 

Multiple Lorentzian profiles provide a good description of 
the powe r spectra in low-mass accretin g neutron star bi naries 
(see e.g. Ivan Straaten et al.1 '2002'. '2003^, 'Reig et al. "2004) and 
black-hole systems (see e.g. Nowak 2000; Belloni et al. 2002; 
iPottschmidt et al.ll2003h . To facilitate the comparison between 
the noise components and timing rapid variability of all types 
of X-ray binaries we have followed this approach here too. The 
peaks from the pulsations were fitted to Lorentzian functions 
with fixed frequency (at the expected value according to the 
spin period and its harmonics) and width (=0.001 Hz w l/T). 
Tables [3] and |4|give the results of the power spectral analysis. 
Figures |8}{9| show the power spectra and noise components of 
various spectral regions. Any given power spectrum can be fitted 
with the sum of two to four broad Lorentzian profiles (excluding 
the peaks), although we identify seven different types of noise 
components. 



4.4. Energy spectra 

The main objective of this work is to investigate the aperiodic 
variability of high-mass X-ray pulsars in correlation with the 
spectral states as defined in CD/HIDs, in the same way as it has 
been extensively done in low-mass X-ray binaries. A detailed 
spectral analysis is out of the scope of this work and will be per- 
formed in a future paper. The reader is referred t o IWilson et al.l 
(2008) (EXO 2030-H375). iMowlavi et alJ (l2006l) (¥ 0332+53). 
Nakajima_eLalJ(|2006|) (4U 1 1 5 +63) and|NaikeLay (l2006|) (KS 
1947+300) for a detailed analysis of the energy spectra during 
type II outbursts. However, for the sake of clarity in the discus- 
sion, it is illustrative to know the X-ray luminosity of the diff'er- 
ent spectral regions. Thus we obtained a time-averaged spectrum 
for each spectral region. The 3-30 keV spectrum was fitted with 
an absorbed power-law model modified at high energies by an 
exponential cutoff. In addition, a Gaussian profile at ~ 6.4 keV, 
that accounts for iron emission line, was required in most cases. 
The X-ray flux and luminosity in the 3-30 keV energy range is 
given in Table |2] The distance to the sources is given in Table [T] 



5. Results 

5.1. Outburst profiles 

The total duration of the 4U 01 15+63 outburst was ~55 days (~ 
MJD 53245-53300) but the rise was shorter (-11 days) than the 
decay (~27 days). The brightest phase of the outburst displays 
a multi-peak profile and lasted for at least 15 days. The source 
spent about one third of the time in this bright state. The maxi- 
mum flux was achieved on ~MJD53264 and amounted to ~280 
mCrab. The PC A observations of 4U 01 15+63 began at the mid 
point of the rising part of the outburst on September 6, 2004 
(MJD 53254) and finished on October 25, 2004 (MJD 53304). 

The ASM light curves of KS 1947+300 and EXO 2030+375 
show very similar characteristics, despite the difference in over- 
all luminosity (the EXO 2030+375 outburst was a factor of 4 
brighter). Not only the outburst profile and duration (~165 days 
and ~155 days, respectively) are similar but also the pre- and 
post-outburst X-ray activity. In these two sources the main out- 



10 



X 

> 



Reig et al.: Rapid variability of Be/X-ray binaries 

KS 1947+300 



0.1 



CO 

E 



I 0.01 h- 



0.001 Ef 



0.0001 






0.01 0.1 1 10 100 



0.01 0.1 1 10 100 



0.01 0.1 1 10 100 



CO 

E 

^ 0.01 



> 



0.001 



0.0001 




0.01 0.1 1 10 100 



0.01 0.1 1 10 100 
Frequency (Hz) 



=1 lllllllj 1 lllllllj 1 lllllllj 1 lllllllj 1 lllllllj ^ 


: 0.1 


|t 


decay4 - 
HB ^ 


— 




+ 


t| = 


\f 






i 






•' 


'1 z 


I / 
; 


II IK 
ii 1 1 ih 

iiiii'.ih 


\ 


* 


^ 




= 


=r 1 1 1 1 


II 


1 h 


\ 1 iiiliiiil-d 



0.01 0.1 1 10 100 



Fig. 6. Power spectra of KS 1947+300 and fit functions of various spectral regions. The number at the top left of each panel is the 
X-ray luminosity in units of 10^' erg s '. HB stands for horizontal branch. The different lines indicate the individual Lorentzian 
components of the fit: L/, (dotted), L/ (thick solid black), L„ (thin solid black). The pulse peaks (fundamental and harmonics) are 
represented by thin dashed orange lines. 



burst is preceded by a minor one, which started one orbital pe- 
riod before the onset of the major event. Likewise, the major out- 
burst was followed by a series of smaller outbursts separated by 
a time interval consistent again with the orbital period (40 days 
in KS 1947+300 and 46 days in EXO 2030+375). The outburst 
profiles are quite symmetric with the rise and decay lasting for 
about 60-70 days. The KS 1947+300 ASM Hght curve contains 
a gap during the rise due to its proximity to the Sun at that time. 

V0332+53 showed the most intense outburst. The maximum 
flux recorded in the energy range 1.3-12.1 keV reached ~1.5 
Crab. The outburst began in 2004 November (MJD 53330) and 
reached maximum flux about one month later (~ MJD 53368). 
The total duration of the outburst was ~ 105 days. As in 4U 
01 15+63, the decay was slower than the rise. The difiference ap- 
peared at the end of the decay, where a longer tail can be seen 
(Fig.©. 

There appears to be a relationship between the shape and du- 
ration of the outbursts and the orbital period of the system. KS 
1947+300 and EXO 2030+375 have the wider orbits and show 
the longer and more symmetric outbursts. The outburst profiles 
of 4U 01 15+63 and V0332+53 have positive skewed (elongated 
tail at the right) and last shorter The skewness statistics (the third 
moment about the mean normalised to the standard deviation, 
jl Yu (^i ~ x)^ /(^^) gives the following values for 4U 0115+63, 
V0332+53, KS 1947+300 and EXO 2030+375: 0.24, 0.58, 0. 16 
and -0.035, respectively. In Be/X-ray binaries, the reservoir of 



matter available for accretion onto the neutron star companion 
comes from the circumstellar disk around the Be star's equator 
The duration of the X-ray outbursts is expected to be propor- 
tional to the amount of accreted material. Thus longer outburst 
duration implies larger disks. This result agrees with the cor- 
relation between the orbital period and strength of the Ho- line 
(iReig et al.lll997l;lReig||2007l) . The longer the orbital period, the 
wider the orbit, the larger the disk. This correlation is explained 
by the truncation of the disk by the neutron star, which will be 
more efficient for narrower orbits. lOkazaki & Negueruelal(l2001l) 
estimated that the Be discs in 4U 01 15+63, V0332+53 and EXO 
2030+375 are truncated at the 4:1 resonance radius, while that 
of KS 1947+300 at the 3:1 resonance. I n terms of stellar radii, 
assuming the stellar parameters given in lOkazaki & Negueruelal 
(i2001i) . these resonance radii translate into truncation radii ~ 4.7 
R., ~ 5.6 R„ ~ 6.9 R, and ~ 8 R. for 4U 01 15+63, V0332+53, 
EXO 2030+375 and KS 1947+300, respectively. This numbers 
should be taken with caution since the calculations are strongly 
dependent on the disc viscosity and assume the same disc phys- 
ical parameters in all systems. 



5.2. Colour analysis 

Be/X-ray binaries show two spectral branches in the CD/HID. 
We shall refer to these branches as the low-intensity or horizontal 
branch and the high-intensity or diagonal branch. However, not 
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Fig. 7. Power spectra of EXO 2030+375 and fit functions of various spectral regions. The number at the top left of each panel is 
the X-ray luminosity in units of 10^^ erg s"' . HB stands for horizontal branch. The different lines indicate the individual Lorentzian 
components of the fit: Li, (dotted), L; (thick solid black), L„ (thin solid black), L,,/ (thick solid grey). The pulse peaks (fundamental 
and harmonics) are represented by thin dashed orange lines. 



all sources display the entire pattern of variability. Sometimes 
only one branch is visible in the CD. In KS 1947+300, the scatter 
of the points that populate the low-intensity branch is too large 
to clearly define the branch (see Fig.|3j. In EXO 2030+375, only 
the hardest (in terms of SC) part of this branch is visible. 

Unless the motion in the horizontal branch is faster than ~ 1 
day ', it seems unlikely that the lack of a complete branch in KS 
1947+300 and EXO 2030+375, at the end of the outburst, is due 
to observational gaps. At least in EXO 2030+375, the frequency 
of the PC A observations was ~ 1 day '. More likely, the lack 
of a clear horizontal branch in the CD of these two sources can 
be attributed to the fact that they did not return to quiescence. 
Both KS 1947+300 and EXO 2030+375 showed smaller out- 
bursts (type I) after the major outburst (type II) analysed in this 
work. The X-ray emission level did not return to zero in these 
two systems (i.e. below the sensitivity threshold of the PCA) as 
can be seen in Fig. |2] The average count rate of the last observa- 
tions was 30 c s' PCU' in EXO 2030+375 and 5 c s' PCU' 
in KS 1947+300. In contrast, 4U 0115+63 seems to have gone 
under the detectability limit of the instruments. It is reasonable 
to think that if the intensity had gone lower, then the horizontal 
branch in KS 1947+300 and EXO 2030+375 would have shown 
up more clearly. Note that V0332+53 also went back to ~ c 
s ' and also displays two distinct branches in the CD. 

In LMXBs, the patterns observed in a CD are often also 
recognizable in the corresponding HID. Sometimes, however. 



spectral branches are more distinct in one of the two diagrams. 
Whether CD or HID presents the cleanest pattern depends on 
source and on the quality of the data. Figure |4]displays the HID 
of the four X-ray pulsars analysed in this work. In contrast to 
the CD, the two branches are clearly seen for all four sources 
when the HC is replaced by the intensity. The junction of the 
two branches occur when the count rate is at ~10% of the peak 
value (-30% in 4U 0115+63). The second column in Table |2] 
indicates the branch in which the source lay during those par- 
ticular observations. For practical reasons, the name "horizontal 
branch" is also used if intensity is plotted on a logarithmic scale 
in the HID, even though in this case this branch does not appear 
horizontal. 

The low-intensity (horizontal) branch is populated by points 
with a low count rate. Since the PCA monitoring began when 
the outbursts were already in progress, the low-intensity branch 
contains points from the end of the outbursts. The only source for 
which data prior to the onset of the outburst are available is KS 
1947+300. As can be seen in Fig.H] the data points (open black 
circles) that corresponds to the very beginning of the outburst 
also lie in the low-intensity branch. Thus we conclude that the 
horizontal branch represents the start and end point of the source 
in its journey through the CD/HID. In terms of X-ray luminosity, 
this branch appears when Lx $ 10^^ erg s~'. 

The evolution of the X-ray colours as the outburst 
evolves, particularly the hard colour (HC), clearly distinguishes 
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Fig. 8. Power spectra of 4U 01 15+63 and fit functions of various spectral regions. The number at the top left of each panel is the 
X-ray luminosity in units of 10^' erg s '. HB stands for horizontal branch. The different lines indicate the individual Lorentzian 
components of the fit: Lb (dotted), L/ (thick solid black), Lif (thick dashed), L„ (thin solid black). The pulse peaks (fundamental 
and harmonics) are represented by thin dashed orange lines. In "peak" and "decayl" the Lb component is affected by the 3 mHz 
QPO, LvLF (see text and Fig.|5l). 



V0332+53 from the other three systems. Unlike, KS 1947+300, 
4U 0115+63 and EXO 2030+375, both the SC and HC in 
V0332+53 decrease as the count rate increases (see Fig|2]i. This 
different behaviour of the HC implies that the motion in the diag- 
onal branch, as the count rate increases, is from the bottom right 
part of the CD to top left part in 4U 01 15+63, KS 1947+300 and 
EXO 2030+375 and from bottom left to top right in V0332+53. 
Another important difference is that the low-intensity branch 
is harder than the high-intensity branch in V0332+53, while it 
is softer in EXO 2030+375, KS 1947+300 and 4U 0115+63 
(Fig. [3j. In this latter case, the horizontal branch corresponds 
to a low/soft (in terms of the hard colour) state that it is not seen 
in low-mass X-ray binaries or black-hole systems. 

After the peak, the source returns following the reverse track. 
The sources spend most of the outburst in the diagonal branch. 
The sources stay in this branch about 70%-80% of the dura- 
tion of the outburst. As can be inferred from Table |2l typical 
timescales are months in the diagonal branch and weeks in the 
horizontal branch. 

4U 0115+63 exhibits hysteresis, that is, the X-ray colours 
adopt different values depending on whether the source is on 
the rise or the decay of the outburst, despite that the count rate 
is similar This effect is also present in V0332+53 and EXO 
2030+375, although not as significant as in 4U 0115+63. It is 
absent in KS 1947+300 (see Fig. gj. In EXO 2030+375 and 



V0332+53 both X-ray colours are larger during the rise than 
during the decay, that is, the sources present a harder spectrum 
during the rise. In 4U 01 15+63, the changes are more complex. 
During the rise SC is larger and HC is smaller than during the 
decay. Although the larger error bars in KS 1947+300 might pre- 
vent us from detecting hysteresis in this source, an effect with 
the amplitude seen in 4U 0115+63 should also be seen in the 
SC-intensity diagram of KS 1947+300. 

Interestingly, cyclotron resonance scattering features have 
been observed for all sources except KS 1947+300. In partic- 
ular, 4U 0115+63 s hows up to five cyclotron Une harmonics 
dCoburn et alj|2004b . However, while the ener gy of the funda 



mental line is found at 1 1 keV in 4U 01 15+63 (Santangel oet al 
1999.) . it appears at 27 keV in V0332+53 (.Pottschmidt et al 



20051). The cyclotron line in EXO 2030-1 -375, at 11 keV, i s 
less significant and lacks harmonic content dWilson et alj|2008h . 
According to the definition of the X-ray colours given in 
Sect. 14.21 the HC in 4U 01 15+63 is expected to be more strongly 
affected by the presence of the cyclotron line and its harmonics, 
while in V0332+53 and EXO 2030+375 its contribution to the 
variability in HC would be insignificant. Moreover, the energy 
of the fundamental cyclotron line in 4U 01 15+63 has been seen 
to increase from ~1 1 kev to 16 k eV as the X-ray lum inosity de- 
creased below 5 X 10^^ erg s ' dNakaiima et al.l 12006) . The re- 
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Fig. 9. Power spectra of V0332+53 and fit functions of various spectral regions. The number at the top left of each panel is the 
X-ray luminosity in units of 10^^ erg s '. HB stands for horizontal branch. The different lines indicate the individual Lorentzian 
components of the fit: Lb (dotted), Li (thick solid black), Llf (thick dashed), L, (dashed-dotted), L„ (thin solid black). The pulse 
peaks (fundamental and harmonics) are represented by thin dashed orange lines. 



lationship between hysteresis, as seen in the HID, and cyclotron 
lines needs to be verified through detailed spectral analysis. 



5.3. Aperiodic variability 

The power spectra of Be/X-ray binaries is dominated by the 
peaks of the pulse component and its harmonics. In order to 
study the aperiodic variability (band-limited noise and possi- 
ble QPOs), the pulse noise should be removed. However, this 
is not an easy task. Some authors substituted the points contam- 
inated by the pulses and their harmonics by a spline fit to the 
remaining points or to a power law conne cting the average of 
a few points before and after the pulse gap dBelloni & Hasingen 
11990 ). This procedure does not give satisfactory results if the 
peaks are too broad. O thers sim ply remove the points from the 
power spectra JReig et al...2006.) . lAngelini et al.l dl989t) removed 
the pulse noise after fitting it with a model consisting of the first 
10 harmonics of the coherent pulsar signal and a power law. 
In this work we have opted to maintain the peaks and fit them 
with Lorentzians. In addition to the distortion introduced by the 
pulse noise, the power spectra continuum of Be/X-ray binaries 
contain significant substructures only visible in certain states. 
Sometime the continuum shows bumps and wiggles, making it 
difficult to obtain good fit s. This difficulty was also found by 
iBelloni & Hasingen (Il990l) . who could not fit the power spectra 
of a number of high-mass X-ray binaries with simple analytical 



laws, such as power law and cutoff, as in low-mass X-ray bina- 
ries. 

In general, the power spectra of Be/X-ray binaries is char- 
acterised by band-limited noise (BLN), that is, a steeping of the 
power density toward high frequency and a flattening toward low 
frequency. Although a double broken power law may give good 
fits in some cases (KS 1947+300) it cannot fit all sources and all 
states (EXO 2030-1-375). In addition, all sources present peaked 
noise in certain spectral states. 

Despite this complexity, we find that the power spectra of 
Be/X-ray binaries can be fitted acceptably with the sum of a 
small number of Lorentzian functions. Ignoring the peaks of the 
pulse component, 2-A Lorentzian profiles are needed to obtain 
acceptable fits. However, we identify 7 different types of noise 
components. 4 of these components are band-limited (BLN) and 
3 QPO. The results of the power spectral fitting are given in 
Tables |3}{5] Figures [8]-|9] show the power spectra of the four X- 
ray pulsars for various regions. Different line types indicate dif- 
ferent noise components as explained below. 

Lh (dotted line) is a zero-centred Lorentzian that describes 
the noise below 0.01 Hz. Due to the very low characteristic fre- 
quency of this component (often peaks outside the available fre- 
quency range) it is not always well constrained. In some states of 
4U 01 15+63, Lb is modified by the appearance of a slow QPO. 
Lb disappears in the horizontal branch of all sources except in 
V0332+53. 
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Fig. 10. Correlation between the characteristic frequencies of the 
noise components. Be/X exhibit lower characteristic frequencies 
than LMXBs. 



Li (thick solid line) is also a broad component with Q < 0.4 
or centred at zero. Its characteristic frequency shifts to low fre- 
quencies when the count rate is low. This component dominates 
most of the power spectrum, covering the range from ~0.1-10 
Hz. In some cases (V0332+53), it is the only component fitting 
the high-frequency end of the power spectrum. 

In 4U 0115-1-63, KS 1947-t-300 and EXO 2030+375 one 
or two more Lorentzians are needed to account for the high- 
frequency noise. We shall call these component L„ (black thin 
solid line) and L,,- (grey thick solid line), with characteristic fre- 
quencies v„ » 1 - 3 Hz, v„' ^ 4 Hz, respectively and Q ~ 0- 0.5. 
These components only appear when the count rate is high, i.e., 
near and at the peak the outburst. 

In general, L/, accounts for the noise short-ward of the main 
pulse peak, while L/ and L„/Lj,' fit the power spectrum above 
the spin frequency. Together they account for the band-limited 
noise. 

The faster pulsars (V0332+53 and 4U 0115+63) display 
peaked noise components, which we shall refer to as QPO noise, 
even though the Q value is not always larger than 2. Table |5] 
gives the best-fit parameters of these components. Llf (dashed 
line in Fig [8] and |9]l is common to the two sources and has a 
characteristic frequency vlf ~ 0.05 Hz and 0.5 > 2lf ~ 2. It be- 
comes stronger as the count rate decreases. Note that although 
the centroid frequency of this Lorentzian component increases 
as the count rate decreases, its width decreases. The characteris- 
tic frequency (i.e. maximum frequency) remains fairly constant 
at the value indicated above (~0.05 Hz). V0332+53 shows, in 
addition, another QPO (that we shall call Lj, dotted-dashed line) 
whose centroid frequency coincides with the frequency of the 
fundamental peak (0.23 Hz) of the pulse period and remains 
constant throughout the outburst. The strength of L, increases 
with count rate, reaching maximum during the peak of the out- 
burst. However, it disappears in low-intensity states, when Lx 
goes below ~ 10^^ erg s"'. 4U 0115+63 also shows another 
QPO at ~ 3 mHz. This QPO is even visible in the light curve 
(Fig. [5]l and has also been reported in previou s outbursts of the 
source dHeindl et alJll999l: ICoburn et al.ll2004 . We have called 
this component LniQPO- 



6. Discussion 

In this section, we compare the correlated spectral and timing 
variability of high-mass X-ray binaries (HMXBs) and low-mass 
X-ray binaries (LMXBs). We comment on the similarities and 
differences of these two types of X-ray binaries. 



6.1. Comparison with low-mass X-ray binaries 

The Be/X-ray binaries that display a wide dynamic range in 
intensity (Type II outbursts) trace two-branch patterns in their 
CD/HID. As in low-mass systems, the sources do not jump 
through the diagram but moves smoothly, following the pattern. 
The horizontal branch (HB) corresponds to a low-intensity state 
and shows the larger fractional nns, similar to the the island state 
in atolls and horizontal branch in Z sources. This branch repre- 
sents the start and end states of the source through the outburst. 
In the high-intensity diagonal branch (DB), the noise compo- 
nents display higher characteristic frequencies and lower nns 
than in the HB. 

The power spectra of LMXBs in the frequency range 
below 100 Hz are c haracte r ised by (Bellon i et al. 200^ 
van Straaten eTall [20021 120031; lOlive et al.1 120031; Ivan der KlS 



2006b : i) a power-law red noise component in the lowest fre- 



quency range, traditionally known as very-low frequency noise 
(VLFN). This component is weak or absent in the island and 
horizontal branches of atoll/Z sources and increases in strength 
as the source moves along the branches; //) strong band-limited 
noise (BLN), consisting of flat-topped noise that breaks at a fre- 
quency Vh ~ 0.01 - 50 Hz. This component can be represented 
by a zero-centred Lorentzian and fits the low-frequency end of 
the power spectrum; ///) peaked noise at frequency ~ 5v/,. When 
this peaked noise appears as a QPO then it is called Lif and its 
characteristic frequency vif. Otherwise, it is normally referred 
as Li,. Sometimes, both components are present, i.e, a narrow 
L^f and a broad L/,; iv) two zero-centred Lorentzians fitting the 
high-frequency part of the power spectrum. They normally ap- 
pear when Vh < 1 Hz and are given the names L/ and L„, respec- 
tively. It is always observed that v/, < vif < v/ < v„. 

In Be/X-ray binaries the power spectra can also be well fitted 
by a relatively small number of Lorentzian components. Some 
of these components can be identified in the power spectra of 
LMXBs. As in LMXBs, the band-limited noise in Be/X-ray bi- 
naries is fitted by three broad Lorentzians: Lb at low frequen- 
cies and Li and L„ at high frequencies. In Fig. [10] we have plot- 
ted the relationship between the characteristic frequencies, i.e. 
Vmax = (vl + (FWHM/2)^y'^, of Lb, Llfm and L/ for the HMXBs 
4U 01 15+63 (fille d squares) and V03 3 2+53 (open squares) and 
for th e LMXBs (JBelloni et al.l 120021; Ivan Straaten et all l2002l 
120031; iReig et al.ll2'004l) 4U 1728-34 (black filled circles), 4U 
1608-52 (cyan filled triangles), GS 1826-24 (orange filled cir- 
cles), IE 1724-3045 (open circles), 4U 0614+09 (diamonds) and 
Aql X-1 (stars). For Aql X-1, L;, was plotted. They are known 
as the WK relation (after lWiinands & van der Klis ( 1999a), up- 
per panel) and the PBK relation (after jPsaltis et al., (il999i) . lower 
panel). The fact that the data points corresponding to the Be/X- 
ray binaries fall on the global correlations found for LMXBs, if 
extrapolated at low frequencies, supports the identification of the 
noise components. 

There are some important differences in the rapid aperiodic 
noise components between HMXBs and LMXBs, namely: /) the 
lower characteristic frequencies, //) the lack of a clear correlation 
of the noise component parameters with mass accretion rate and 
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Hi) the noise component, Lj, associated with the peaks of the 
periodic modulation. 

In accretion-powered X-ray pulsars, QPOs fall in the milliHz 
range, that is few order of magnitudes lower than in LMXBs. 
There are about 12 high-mass X-ray pulsars that show QPOs 
in their power spec tra. The frequencies of t hese QPOs range 
from 1-400 mHz ('Shirakawa&Tai' '2002'; 'Inam etalj l2004t 
iMukheriee eTal. 2006) . Recently, Kaur et al. (2007) reported the 
discovery of a 1.3 Hz QPO in XTE JOll 1.2-7317, a Be/X-ray 
binary located in the SMC. Many of these QPOs have Q values 
lower than 2, hence strictly speaking, they should be referred as 
peaked noise (conventionally, peaked noise is considered as a 
QPO when the quality factor Q - v/FWHM > 2). However, 
since these features (also L^f and L., in this work) represent a 
substantial concentration of power in a limited frequency range 
they are referred to as QPOs. 

If it is assumed that QPOs are formed as a result of Keplerian 
motion of inhomogeneities in an accretion disk then the longer 
implied timescales of the QPOs in HMXB is somehow ex- 
pected as the inner radius of the accretion disk must be larger 
than the magnetospheric radius (r,„ ~ 10** cm). 10-400 mHz 
QPOs imply disk radius of the order of 10** - 10^ cm {Rj ~ 
(GM I AtP-v^pq)^ ^^). The origin of the mHz QPO is, however, un- 
clear. N either the Keplerian frequency model ([van der KUs et alj 
I1987 1J) nor the beat frequency model JAlpar & Shaham"l985') 
can explain the origin of these QPO in all systems. QPOs be- 
low the neutron star spin frequency cannot be explained by the 
Keplerian frequency model because if the Keplerian frequency 
at the magnetosphere is less than the spin frequency then th e 
propeller mechanism would inhibit accretion (S tella et aljl986h . 
This is the case of 4U 0115+63 and V0332-H53. Likewise, the 
beat frequency model predicts a centrifug al inhibition t hreshold 
that it is at variance with the observations (Finger"1998'). A third 
model, the magnetic disk precession model ( Shirakaw a & Lail 
l2002h . attributes the mHz QPO in X-ray pulsars to warp- 
ing/precession modes induced by magnetic torques near the in- 
ner edge of the accretion disk. The applicabihty of th is model to 
all sources is not as yet certain ( Mukheriee et al.ll20() 6). 

Unlike LMXBs where the power spectral parameters and 
strength of the various noise components (anti)correlate with the 
mass accretion rate, M, HMXB noise appears to depend weakly 
on M (we assume that the accretion rate increases/decreases dur- 
ing the rise/decay of the outburst), at least at high and intermedi- 
ate accretion rates. Neither the characteristic frequencies nor the 
rms of the Lorentzians that account for the band-limited noise 
show a clear relationship with X-ray flux in the diagonal branch, 
despite the change of about 1 .5 order of magnitude in luminosity 
(see Tables[3H5]). The lack of such correlati ons were also noticed 
in the 1985 outburst of E XO 2030+375 dAngelini et al. 1989r: 
iBelloni & Hasingedll990h . 

Ls is exclusive of accreting X-ray pulsars and it is associ- 
ated with the pulse noise. Its central frequency coincides with 
the frequency of the fundamental peak of the pulse period and 
suggests a strong coupling between the periodic and aperiodic 
noise components. This coupling implies that the instabilities in 
the accretion flow that give rise to the aperiodic variability must 
travel all th e way down to the neut ron sta r surface. The read er is 
refen-ed to lLazzati & Stellal (Il997h and lBurderietal.1 (1 1997b for 
a detailed study of the coupling between the periodic and ape- 
riodic variability in X-ray pulsars and to iReig et all (l2006h for 
the particular case of V0332+53. L^ is absent in the horizontal 
branch. In the case of V0332+53, the disappearance of the pulse 
peaks in the horizontal branch is accompanied by the disappear- 
ance of Ls, which supports the idea of a common physical origin. 
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Fig. 11. Representative power spectra of high- and low-intensity 
states. Note the different shape of the noise at low frequencies. 
Hat-topped noise is seen in the HB of KS 1947+300 (and EXO 
2030+375 and 4U 01 15+63) and p ower-law noi se in V0332+53. 
Confront this figure with Fig. 4 of Iva n der Khsl (11994) and note 
the similarities between KS 1947+300 and atoll sources and be- 
tween V0332+53 and Z sources. The "peak" power spectrum 
was shifted by x 100 for easy display. 



6.2. Are there atoll and Z sources in HMXBs? 

The correlated timing and spectral analysis defines two different 
types of Be/X-ray binaries. The differences are not only seen in 
the shape of the CD but also in the number and properties of the 
noise components. 

The pattern traced by V0332+53 in the CD reminds that of 
the Z-sources, with the diagonal branch being the analogue to 
the normal branch and the HMXB horizontal branch branch be- 
ing the counterpart of the LMXB horizontal branch. The flaring 
branch would be missing in V0332+53. As in Z sources, the 
hardest spectrum corresponds to a lower count rate state. Also, 
the fractional rms is higher and the characteristic frequencies are 
lower in the horizontal branch than in the diagonal branch. 

In 4U 01 15+63, KS 1947+300 and EXO 2030+375, the HB 
also corresponds to a low-intensity state and also displays the 
larger rms of the noise components. However, the HB in these 
three sources is softer than the diagonal branch. Likewise, the 
overall shape and motion along the branches diff'ers from that 
seen in V0332+53 and LMXBs atoll sources. The HB in 4U 
0115+63, KS 1947+300 and EXO 2030+375 corresponds to a 
low-intensity spectrally soft state not seen in other types of X- 
ray binaries. 

Further similarities between V0332+53 and Z sources and 
4U 0115+63, KS 1947+300 and EXO 2030+375 and atolls 
sources are found in the shape of the very-low frequency noise 
(<0.1 Hz). Fig, im shows the power spectra of the HB and peak 
of the outburst for V0332+53 and K S 1947+300 This figure 
should be compared to, e.g.. Fig. 4 in Ivan der Klisl ( 11994 ). The 
noise below -0.1 Hz in the HB is flat-topped in KS 1947+300 
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(and in EXO 2030+375 and 4U 01 15+63), as in the island state 
of atoll sources, and power-law like in V0332+53, as in the HB 
of Z sources. 

Initially, the differences between LMXB Z and atoll sources 
were ascribed to a higher magnetic-field strength and higher 
accretion rates in Z sources. However, the detection in atoll 
sources of QPOs that are reminiscent of HBO an d NBO 
dFord & van der Klisl 119981: IWiinands & van der Khsl [l999b) 
and the existence of sources that display t he two types of be- 
haviour, Z and atoll dHoman et al.l l2007alfbr) . question the va- 
lidity of the proposed higher magnetic field in the Z sources. 
Nevertheless, it is worth noting that V0332+53 displays the most 
intense outburst and contains the stronger magnetic field of the 
four sources. The peak luminosity in V0332+53 is about a fac- 
tor 2, 3 and 5 that of EXO 2030+375, 4U 0115+63 and KS 
1947+300, respectively. The strength of the magnetic field can 
be estimated from the energy of cyclotron resonant scattering 
features (CRSF) through the relation 



■^cyc 



11.6 



B 



1012G 



(l+zr keV 



where B is the magnetic field and z is the gravitational redshift. 
As mentioned above, KS 1947+300 is the only source for which 
no CRSF has been reported, while for 4U 0115+63 E^y^ - 12 
keV (l^ntangelo et al. 1999), for EXO 2030+375 £cyc = H 
keV (Wilson etal. 2008), and for V0332+53 £cyc = 27 keV 
JPottsch midt et al. 2005). Thus the magnetic field in V0332+53 
is about 2-3 times stronger than in the other three systems. 

It is tempting to attribute the differences between Be/X-ray 
binaries to the source brightness and magnetic field strength. 
The power-law noise at low frequencies in the HB of V0332+53 
could be due to the larger intensity and the presence of Ls noise 
to a higher magnetic field. 4U 0115+63, KS 1947+300 and 
EXO 2030+375 show flat-topped noise and lack L, component. 
However, with only one source of this type no definitive con- 
clusion can be drawn about whether the source brightness and 
magnetic field strength allows the distinction between different 
types of systems among HMXB. Future observations will clarify 
this issue. 



7. Conclusions 

We have performed, for the first time, a systematic analysis of 
the spectral (through the study of hardness ratios) and rapid ape- 
riodic variability of four Be/X-ray pulsars during type 11 out- 
bursts. Our aim was to characterise the rapid aperiodic variabil- 
ity in correlation with the position in the colour-colour diagram 
as it has been done for black-hole systems and low-mass X-ray 
binaries. 

The correlated spectral-timing behaviour in Be/X-ray bina- 
ries shares a number of similarities with low-mass X-ray bina- 
ries, namely, /) existence of spectral branches in the CD/HID. At 
high and intermediate flux the sources move along the diagonal 
branch; at very low count rate the soft colour decreases, while 
the hard colour remains fairly constant defining a horizontal 
branch //) smooth motion, i.e. without jumps, in the CD/HID, ///) 
the low-intensity states are more variable (in terms of fractional 
rms), iv) description of power spectra in terms of a small num- 
ber of Lorentzian components, v) similar BLN namely, broad 
Lorentzians Q ~ Q- 0.3 describing the low-frequency and high- 
frequency noise (Li,, Li, Lu) and peaked noise in between (Lif) 
and vi) flat-topped noise at the lower frequencies in the horizon- 
tal branch that turns into power-law noise in the diagonal branch. 



There are also important differences: /) different patterns 
in the CD/HID. Unlike low-mass X-ray binary and black hole 
systems, the low-intensity state is associated with a soft state 
(except for V0332+53), ii) slower motion along the spectral 
branches, hours to days in LMXBs, weeks to months in HMXBs, 
///) the characteristic time scales implied by the noise compo- 
nents are about one order of magnitude longer in HMXBs (e.g. 
mHz QPOs), iv) Ls noise in HMXBs and v) no apparent cor- 
relation between the power spectral parameters (characteristic 
frequencies, rms) and mass accretion rate. 

Although we have analysed all Be/X-ray binaries that 
showed type 11 outbursts in the period 1996-2007 and had good 
data coverage and sampling, the number of sources is not large 
enough to draw definite conclusions on the underlying physi- 
cal parameters that may explain the differences between the two 
subgroups. Increasing the number statistics is, however, difficult 
given the transient nature of their X-ray emission and unpre- 
dictabihty of the type II outbursts. 

Even though all four sources analysed in this work are X- 
ray pulsars, it is not clear whether the high-mass X-ray binaries 
with supergiant companions show a similar behaviour In these 
systems accretion occurs via the strong stellar wind of the mas- 
sive companion, which contrasts with accretion via a disc, as it 
is believed to occur in Be/X-ray binaries during outburst. The 
problem of this type of studies in supergiant X-ray binaries is 
that they are rather stable in long time scales and do not show a 
wide range in X-ray flux. 

Acknowledgements. This reseai'ch was supported by the European Union Marie 
Curie grant MTKD-CT-2006-039965. This research has made use of NASA's 
Astrophysics Data System Bibhographic Services and of the SIMBAD database, 
operated at the CDS, Strasbourg, France. The ASM hght curve was obtained 
froin the quick-look results provided by the ASM/RXTE team. 



References 

Angelini, L.; Stella, L.; Parmar A. N., 1989, ApJ, 346, 906 

Alpar, M. A. & Shaham, J. 1985, Nature, 316, 239 

Belloni, T, & Hasinger, G. 1990, 230, 103 

Belloni, T, Psaltis, D. & van der Klis, M. 2002, ApJ, 572, 392 

Belloni, T, Homan, J., Casella, P. et al. 2005, A&A, 440, 207 

Bernacca PL., lijima T, Stagni R. 1984, A&A, 132, L8 

Borozdin, K., Gilfanov, M., Sunyaev, R., et al., 1990, Sov. Astron. Lett., 16, 345 

Burderi, L., Robba, N.R., La Barbera, N., Guainazzi, M. 1997, ApJ, 481, 943 

Cainero Arranz, A., Wilson, C. A., Coimell, P, et al. 2005, A&A, 441, 261 

Chakrabarty, D., Koh, T., Bildsten, L., et al. 1995, ApJ, 446, 826 

Cobum, W., Kalemci, E., Kretschmar, P, et al. 2004, ATel 337 

Coe M.J., Longmore A.J., Payne B.J., Hanson C.G. 1987, MNRAS, 226, 455 

Coe, M. J., Payne, B. J., Longmore, A., Hanson, C. G., 1988, MNRAS, 232, 865 

Corbet R.H.D., Charles PA., van der Klis M. 1986, A&A, 162, 1 17 

Finger, M. H., Wilson, R. B., Harmon, B. A. 1996, ApJ, 459, 288 

Finger, M. H., 1998, AdSpR, 22, 1007 

Ford, E. & van der Klis, M. 1998, ApJ, 506, L39 

Galloway, D.K., Morgan, E.H., Levine, A.M., 2004, ApJ, 613, 1164 

Goranskii, V.P, Esipov, V.F., Lyutyi, V.M., & Shugarov, S.Yu., 1991, Sov Astron. 

Lett., 17, 399 
Goranskij, V. & Barsukova, E. 2004, ATel 245. 
Grankin, K.N., Shevchenko, VS., Yakubov, S.D., 1991, Sov. Astron. Lett, 17, 

415 
Hasinger, G., 1987, A&A, 186, 153 
Hasinger, G., van der Klis, M., 1989, A&A, 225, 79 
Hayasaki, K. & Ozakaki, A.T. 2004, MNRAS, 350, 971 
Heindl, W. A., Coburn, W., Gruber, D. E., Pelting, M. R., Rothschild, R. E. et al. 

1999, ApJ, 521, L49 
Heindl, W. A., Rothschild, R. E., Coburn, W. et al., 2004, in X-ray Timing 

2003: Rossi and Beyond. AIP Conference Proceedings, Vol. 714, held 3- 

5 November, 2003 in Cambridge, MA. Edited by Philip Kaaret, Frederick 

K. Lamb, and Jean H. Swank. Melville, NY: American Institute of Physics, 

p.323-330 
Hoinan, J., van der Klis, M. Wijnands, R. Belloni, T, Fender, R. 2007a, ApJ, 

656, 420 



Reig et al.: Rapid variability of Be/X-ray binaries 17 



Homan, J., Wijnands, R., Altamirano, D. Belloni, T. 2007b, ATel 1 165 

Inam, S. ., Baykal, A., Swank, J., Stark, M. J., 2004, ApJ, 616, 463 

Jemigan, J.G., Klein, R.I., Arons, J., 2000, ApJ, 530, 975 

Kaur, R., Paul, B., Raichur, H., Sagar, R., 2007, ApJ, 660, 1409 

Klochkov, D., Horns, D., Santangelo, A., et al. 2007, A&A, 464, L45 

Kriss, G. A., Cioffi, D. R, Canizares, C. R. 1983, ApJ, 272, 439 

Lazzati, D. & Stella, L., 1997, ApJ, 476, 267 

Leahy, D.A., Darbro, W., Eisner R-F. et al., 1983, ApJ, 266, 160 

Levine, A., & Corbet, R., 2000, lAUC No 7524 

Makishima K., Mihara T., Ishida M., et al. 1990, ApJ, 365, L59 

Markert, T.H, Canizares, C.R., Clark, G.W., et al. 1977, ApJ, 218, 801 

Mavromatakis, F., 1994, A&A, 285, 209 

Miyamoto, S., Kimura, K., Kitamoto, S., Dotani, T., Ebisawa, K. 1991, ApJ, 383, 

784 
Mowlavi, N.; Kreykenbohm, 1.; Shaw, S. E. et al. 2006, A&A, 451, 187 
Motch, C, Janot-Pacheco, E., 1987, A&A, 182, L55 

Motch, C, Stella, L., Janot-Pacheco, E., & Mouchet, M. 1991, ApJ, 369, 490 
Mukherjee, U., Bapna, S., Raichur, H., Paul, B., Jaaffrey, S.N.A. 2006, JApA, 

27,25 
Nakajima, M., Mihara, T., Makishima, K., Niko, H., 2006, ApJ, 646, 1 125 
Nowak, M.A., Vaughan, B.A., Wilms, J., Dove, J.B., Begelman, C. 1999, ApJ, 

510, 874 
Nowak, M.A. 2000, MNRAS, 318, 361 

Naik, S., Callanan, P.J., Paul, B., Dotani, T., 2006, ApJ, 647, 1293 
Negueruela I., Roche R, Fabregat J., Coe M.J. 1999, MNRAS, 307, 695 
Negueruela, I., Israel, G.L., Marco, A., Norton, A.J., Speziah, R., 2003, A&A, 

397, 739 
Okazaki, A., Negueruela, 1., 2001, A&A, 377, 161 
Olive, J., Barret, D., Gierlinski, M., 2003, ApJ, 583, 416 
Parmar, A. N., White, N. E., Stella, L., Izzo, C, Ferri, P, 1989, ApJ, 338, 359 
Pottschmidt, K., Wilms, J., Nowak, M. et al. 2003, A&A, 407, 1039 
Pottschmidt, K., Kreykenbohm, I., Wilms, J., et al., 2005, ApJ, 634, L97 
Priedhorsky, W, Hasinger, G, Lewin, W.H.G, 1986, ApJ 306, L91 
Psaltis, D, Belloni, T., van der Klis, M. 1999, 520, 262 
Qu J.L., Zhang S., Song L.M., Falanga, M. 2005, ApJ, 629, L33 
Rappaport, S., Markert, T., Li, F.K., et al. 1977, ApJ, 217, L29 
Rappaport, S., Clark, G. W., Cominsky, L., Li, R, & Joss, P C. 1978, ApJ, 224, 

LI 
Reig, R, Fabregat, J., Coe, M.J., 1997, A&A, 322, 193 
Reig, P, Coe, M. J., 1998, MNRAS, 294, 118 

Reig, P, Stevens, J. B., Coe, M. J., Rabregat, J., 1998, MNRAS, 301, 42 
Reig, P, Coe, M. J., 1999, MNRAS, 302, 700 
Reig, P, van Straaten, S. & van der KUs, M. 2004, ApJ, 602, 918 
Reig, P, Martinez-Niinez, S. Reglero, V., 2006, A&A, 449, 703 
Reig, P, 2007, MNRAS, 377, 867 
Reig, P.,Larionov, V., Negueruela, 1., Arkharov, A. A., Kudryavtseva, N. A., 

2007, A&A, 462, 1081 
Reynolds, A. P, Parmar, A. N., White, N. E., 1993, ApJ, 414, 302 
Santangelo, A., Segreto, A., & Giarrusso, S. 1999, ApJ, 523, L85 
Schulz, N.S., Hasinger, G. Truemper J., 1989, A&A, 225, 48 
Shirakawa, A. & Lai, D. 2002, ApJ, 565, 1134 
Stella L., White N.E., Davelaar J., et al. 1985, ApJ, 288, L45 
Stella, L., White, N. E., Rosner, R. 1986, ApJ, 308, 669 
StoUberg, M.T., Finger, M.H., Wilson, R. B., et al. 1999, ApJ, 512, 313 
Sun, X.J., Li, T.P, Wu, M., Cheng, L.X., 1994, A&A, 289, 127 
Swank, J., Morgan, E., 2000, lAUC No 7532 

Takeshima T., Dotani T., Mitsuda K., Nagase R 1994, ApJ, 436, 871 
Tanaka, Y. 1983, lAUC 3891 

Terrell J., & Priedhorsky W.C. 1984, ApJ, 285, L15 
van der Klis, M., Jansen, R, van Paradijs, J. et al., 1987a, ApJ, 313, L19 
van der Klis, M., Stella, L., White, N., Jansen, R, Parmar, A. N. 1987b, ApJ, 

316,411 
van der KUs, M., 1989, in Timing Neutron stars, NATO ASl ser C vol. 262, ed. 

H.Ogelman,& E.P.J, van den Heuvel, p. 27 
van der Klis, M. 1994, ApJSS, 92, 511 
van der Klis, M. 2006, in Compact Stellar X-ray sources, eds. W.H.G. Lewin & 

M. van der Klis, Cambridge University Press, p. 39 
van Straaten, S., van der KUs, M., di Salvo, T. & Belloni, T. 2002, ApJ, 568, 912 
van Straaten, S., van der Klis, M., & Mendez, M. 2003, ApJ, 596, 1 155 
Wilson, C.A., Ringer M.H., Coe, M.J., Laycock, S., Rabregat, J. 2002, ApJ, 570, 

287 
Wilson, C.A., Ringer, M.H., Camero Arranz, A. 2008, ApJ, 678, 1263 
White, N.E., Marshall, RE., 1984, ApJ, 281, 354 
Whitlock L. 1989, ApJ, 344, 371 
Wijnands,R. & van der Klis, M. 1999a, ApJ, 514, 939 
Wijnands, R. & van der Klis, M. 1999b, ApJ, 522, 965 



